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Abstract 

The aim of this paper is to prove the existence of weak solutions to the equation Ait+it*' = 
which are positive in a domain Q C R^, vanish at the boundary, and have prescribed 
isolated singularities. The exponent p is required to lie in the interval {N/{N — 2), (A'^ + 
2) / (A'^ ~ 2)). We also prove the existence of solutions to the equation Au + = which are 
positive in a domain Q C R" and which are singular along arbitrary smooth fc-dimensional 
submanifolds in the interior of these domains provided p lie in the interval {{n — k)/{n — k — 
2), (n — fc + 2)/(n — fc — 2)). A particular case is when p — (n + 2)/(n — 2), in which case 
solutions correspond to solutions of the singular Yamabe problem. The method used here 
is a mixture of different ingredients used by both authors in their separate constructions of 
solutions to the singular Yamabe problem, along with a new set of scaling techniques. 
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1 Introduction and statements of main results 

In this paper we construct solutions with prescribed singularities for the semilinear elliptic equa- 
tion Au + w^ = and other closely related equations, for a certain range of values of the exponent 
p, in a variety of situations. The solutions will have singularities prescribed along a disjoint union 
of submanifolds of varying dimension. We now describe our results, starting with the simplest 
case, when the prescribed singular set is discrete. 

Suppose is any bounded open set in M.^ , with smooth boundary. Consider the equation 

—Au = in , , 

M = on dn, ^ ' 

A weak solution of (|^) is a function u which solves this equation on all of fl in the sense of 
distributions. In particular, any such solution must belong to U^{Q,). The singular set of a weak 
solution u, sing (u), is the complement in of the set of points where u is continuous, and hence 
smooth. 

Our first result concerns the case when the singular set is finite. We single it out since its 
proof is slightly simpler than the more general case. 

Theorem 1 Let D, be as above, and suppose that E — {xi, . . . , xk} C fl is any finite set of 
points. Suppose also that the exponent p lies in the range 

N N + 2 

<P<T^- (2) 



N -2 N -2 

Then there is a K parameter family of positive weak solutions u of with sing(u) — S. In fact, 
the solution space of this equation is locally a K dimensional real analytic variety. 

It is known that if the exponent p is less than N/{N — 2), then any weak solution of 
must be smooth on all of fl. The existence of solutions of this equation with prescribed isolated 
singularities when p lies in the interval N/{N — 2) < p < po, where po is some value close 
to N/{N — 2) (and in particular, less than {N + 2)/{N — 2)), has already been solved by the 
second author in |l^ and in ||ll|. When p ~ {N + 2)/{N — 2), which is the so-called critical 
exponent, the problem becomes conformally invariant. There is then a loss of compactness and 
the problem consequently becomes much more difficult. Solutions now correspond to metrics of 
constant positive scalar curvature which are complete in a neighbourhood of the singular points. 
It is geometrically more natural in this case to replace the domain by (or in fact, any other 
compact manifold of positive scalar curvature); the operator A then needs to be replaced by the 
conformal Laplacian A — -^^-^Rq, where i?o is the scalar curvature of the background manifold. 
An additional source of difficulties in this case is that the position of the singularities is no 
longer necessarily arbitrary. The general existence result for this geometric problem (when the 
background manifold is the sphere) was obtained by R. Schoen [Q . The recent work of the first 
author, along with D. Pollack and K. Uhlenbeck examines the moduli space of solutions for 
this problem. 

In any event, Theorem |l] extends the result of |l^ and jl^ to the full range of subcritical 
exponents. 

It is also possible to prove existence of solutions of (|l]) which are singular along submanifolds 
of higher dimension. Now let E = U^j^S^, where each C is a smooth submanifold without 
boundary of dimension fc, > 0. We shall suppose that il C M", since we reserve the symbol N 
(or more properly Ni = n — ki) for the dimension of the normal space to each E^. 
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Theorem 2 Suppose that the exponent p satisfies 



n 

<p < — 



hi 2 
or in other words 

2p + 2 2p 
n — < Ki < n — 



p — 1 P ^ 1 

for i = 1, . . . , K . Then there is a positive weak solution u of with sing{u) ~ E. Provided at 
least one of the ki > 0, there is an infinite dimensional space of solutions for this problem. 

Again, the existence of solutions for this equation with higher dimensional singular sets of 
dimension k has been solved by Y. Rebai jisj when the exponent p lies in the interval (n — 
k)/{n — k — 2) < p < Pk, where pk is close to the lower cndpoint of this interval, and in particular, 
less than (n — fc + 2)/(n — fc — 2). We shall only prove this result when p > [n — ki)/{n — ki — 2) 
since it would complicate the notation to cover the borderline case where p attains the lower 
limit of this interval of values. However, the arguments below could be extended without undue 
difhculty to cover this value as well. 

A special case of Theorem ^ is the singular Yamabe problem. This occurs when the exponent 
p attains the critical value (n + 2)/(n — 2). In this case, the dimensions of the Si are allowed to lie 
in the range < ki < (n — 2)/2. As discussed above, it is more natural now to replace 17 by the 
round sphere (S",go), or indeed by any compact manifold (Af, go) with (constant) nonnegative 
scalar curvature, and to replace the Laplacian by the conformal Laplacian Lq = A — -R(go)- 
The relevant equation now is 

- Lqu = (3) 

on M . At least when the background manifold is the sphere, it follows from work of R. Schoen 
and S. T. Yau [|l5j that solutions to this singular Yamabe problem (for which the corresponding 
conformally related metrics g = m'*/^"~^^(7o are complete and have bounded Ricci curvature) exist 
only if the dimension of the singular set is less than or equal to {n — 2)/2. Thus, Theorem || 
implies 

Theorem 3 Let {M, go) be any compact manifold with constant nonnegative scalar curvature. Let 
Yj <Z M be any finite disjoint union of smooth suhmanifolds of dimensions ki with < ki < . 
Then there is an infinite dimensional family of complete metrics on M \ Ti with constant positive 
scalar curvature. 

Together with the result of jl^, which treats the case ki = ^^^^ when n is even, this theorem 
settles the question of existence of solutions to the singular Yamabe problem (with constant 
positive scalar curvature) whenever E is a finite disjoint union of smooth submanifolds with 
dimensions greater than zero, but less than or equal to {n ~ 2)/2. 

Finally, we may apply the arguments of ||^ to slightly refine Theorem ^ 

Theorem 4 Let S be any finite disjoint union of C^'" submanifolds in fl of dimensions ki sat- 
isfying the restrictions above. Then has an infinite dimensional family of solutions. 

The basic idea of Theorems |^ and is that solutions of (|l]) with positive dimensional singular 
sets may be obtained as perturbations of solutions to the problem on the fibres of the normal 
bundle of S. The exponent ^^^1 is subcritical for the induced problem on these fibres, which are 
of lower dimension than the ambient space. In particular, the problem with critical exponent 
may be reduced to a subcritical one. 

The paper is organized as follows. First we analyze the asymptotics of rotationally invariant 
solutions of the problem on \ {0}. These are then used to construct approximate solutions 
for each of the problems. At this point we then give a more detailed outline of the strategy 
of the proof. The main task is to analyze the linearization of the nonlinear operator around 
these approximate solutions, and to prove that these linearizations are surjective on appropriate 
function spaces. This occupies the bulk of the paper. After this is accomplished, the exact 
solutions are obtained by a fixed point argument. The deformation spaces for the solutions of 
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these operators, and the arguments necessary to replace $7 by a manifold M are discussed in the 
last section. 

The authors are grateful to Dan Pollack for having carefully read the manuscript and making 
a number of suggestions to improve the exposition. 

2 Singular radial solutions on \ {0} 

In this section, we recall some well known facts which appear for example in [||. 

Proposition 1 For any exponent p G { ' jv^2 )^ there exists a one parameter family of weak 
solutions «£, e > 0, for the equation 

-^Au^uP in M^, (4) 
such that the are radial, singular only at the origin, and satisfy the following properties : 

• > for < r < oo. 

• Me can be written as 

Ue{x) = |a:|"~t;i(-log(|x|/e)), 
where the function vi is bounded independently of e. 

• limt^+oo Vi{t) = Voo > 0, where 

and in particular, is independent of e. 

• limt^_oo e^*^^^ p-^\i{t) < +Qo; in particular, for \x\ large 

u,{x) = c(e^-F^ + 0{\x\-^)\x\^-'^ , 
where the constant c is independent of e. 

• Finally, WviWlZ} < ^^vP^\ 

Proof : First define a new independent variable t = — log \x\ and set 

u{x) = \x\~~v{—\og{\x\)). (6) 
Then the new function v{t) satisfies 

duV^(N-P±^)dtV-^(N-^)v + vP = 0. (7) 

V J3-1/ p-iy J5-1/ 

Now look at the phase-plane portrait for this equation in the {v,vt) plane. The two equilibrium 
points are (0,0) and {voo,0), where Voo is defined in (||); the first of these is a saddle point and 
the second a stable equilibrium. There is a single orbit issuing from (0,0) and tending to (woo,0) 
as t ^ oo. Let vi(t) be one of the functions corresponding to this orbit; it is determined only up 
to the choice of its initial Cauchy data {vi{0),v[{0)). Now let Ve{t) = vi{t + loge). As e varies 
over (0, oo) we obtain all solutions of (Q) which converge to as i ^ — oo and to Voo as t ^ oo. 

We now obtain an upper bound on the function vi{t). It is easy to see that the trajectory 
{vi{t) , v'i{t)) is contained within the homoclinic orbit of the Hamiltonian system 

dttw - (n - w + wP = (8) 

p-1 \ p-lj 
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which tends to (0,0) as t tends both to +cxd and —oo. Let {'Wi{t),'w[{t)) parametrize this orbit. 
Then we conclude that 

sup vi < sup Wi. 
The conservation of Hamiltonian energy for (|^) now shows that 

1 ,„ . >2 _ _P (,,_ JP A .,,2 j_ J_...r*i 



_(a.,)^__^(^iv--f^j.; + _.r=o. (9) 

wi attains its supremum when w'l = 0, so we obtain the upper bound 

i;P-i(t) < ^-i-^w^-i for all t e M. 
This ends the proof of the Proposition. □ 

Remark 1 For any e > 0, whenever u{x) is a solution of (Qj, e~ u{e~^ x) is another solution 
of this same equation. This dilation invariance corresponds exactly to the translation invariance 
of the equation ^ for v. We let Ui{x) be the solution corresponding to Vi and 

Uf{x) = e p-iui(e^ x) 
Later we will also need the following : 

Remark 2 By this resettling we can always assume that, for any given constant a > 0, vi may 
be chosen so that 

supvi{t)<a. (10) 
t>o 

The final remark relates the construction of approximate solutions here to those used in |^ : 
Remark 3 The stable point (woo,0) corresponds to the singular solution 

2 

Uoix) EE Voo\x\^~- 

This solution is invariant by the scaling described above. The fact that ui is not dilation invariant 
is important in the construction of approximate solutions. 

3 Function spaces 

In this section we define the weighted Holder spaces C^'"{fl \ S) appropriate for this problem. 
Roughly speaking, the functions in these spaces are products of powers of the distance to S with 
functions whose Holder norms are invariant with respect to scaling by dilations from any point 
on S. 

We shall use local Fermi coordinates around each component of S to define these spaces. 
When is a point, these are simply (geodesic) polar coordinates (r, 9) around that point. When 
Si is higher dimensional, let T}^"^ be the tubular neighbourhood of radius cr around S^. It is well 
known that T^'*"* is a disk bundle over S^, and is diffeomorphic to the disk bundle of radius a in the 
normal bundle NT,i. Using the metric, this diffeomorphism is canonical. The Fermi coordinates in 
this tubular neighbourhood will be constructed as coordinates in the normal bundle, transferred 
to rj*^ via this fixed diffeomorphism. Here r is the distance to E^, which is well defined in T^^*"* 
and smooth away from E^ provided a is small enough, y is a local coordinate system on E^, and 
9 is the angular variable on the sphere in each normal space Ny'Si. Let Bn^^ denote the ball of 
radius a in Ny'Si. We shall let x denote the rectangular coordinate in these normal spaces, so 
that r = \x\ and 9 = x/\x\. 

Let us also fix a function p > in C°° (fl \ E) with p equal to the polar distance r in each T^*-*-* . 
Let w be a function in this tubular neighbourhood, and define 

I I I (r + r)°'\w(z) - w(z)\ 

-U) - ^ sup \w\+ sup ^ J I \ J WI 

.Jh ._rjT^^Ar~r\- + \y~y\- + {r + fr\9-9\<- 

Here z, z are any two points in and (r, y, 0), (f, y, 9) are their Fermi coordinates. 
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Definition 1 The space Cq'°'{^ \ S) is defined to be the set of all w G C'^'"(f2 \ E) for which the 
norm 

K k 
i=l j=0 

is finite. Here ^aji = ^'^^'^j for any 7 G M, 

C:^'"(0 \i:) = {w = p^w:we Co'"(n \ s)}. 

Thus functions in Cy"(f7 \ S) are allowed to blow up like near the E^. Note that when is 
positive dimensional, functions in C^'" may be differentiated in the 'tangential' direction only at 
the expense of giving up a power of p. Equivalently, their derivatives with respect to up to fc-fold 

products of the vector fields rdr, rdy, de blow up no faster than p'^. 
We collect a few essentially trivial remarks about these spaces : 

Lemma 1 1- Ifl& (-iV,0) then Cl^'"{n \ E) c LP{n) for p & (1, -N/^). 

2. For w e C^'" and v e C^f, then wv e C^!^^, and also 

,a,7 • 

3. Given w € C^'" with w>0 and if p > k + I, then £ C^^" and 

\\wP\\k,a,in < CpMll^c.^^, 
for some constant Cp > only depending on p. 

4. If we C^+i'° and \Vw\ e C^:°i, then for any a e (0, 1), w G C^'«. 

We shall also need a slightly better estimate for the norm of w^. 

Lemma 2 For any 77 > 0, 7 G M, there exists some > 0, depending only on 7 and rj, such that 
if \\w\\o,a,j < 6 then 

||wP||o,a,7-2 < ??||w||o,a,7- 

Proof : The proof relies on the simple fact that 

I Kilo ,a,P7 ^ erf ||w||o,a,7- 
Since 7 > we see that P7 > 7 — 2, and the result follows immediately. □ 

Definition 2 

C^'^(17 \ E) = {u; G C:5'"(17 \ E) : w = on d^}. 

Finally, we shall also need to use Holder spaces on \ {0}, respectively M" \ R*^, which 
have different decay properties at the origin, respectively at R'', than at 00. Now we simply use 
the global polar coordinates (r, 9) in \ {0}, and cylindrical coordinates (r, 9, y) in R" \ R*^ = 
R^ \ {0} X R'=. 

Definition 3 For any 7,7' G R, the space C^'y(R-'^ \ {0}) consists of all functions w for which 
the norm 

||w'||fe,a,7,7' = sup + ^Up ||w||fe,a,^/ 

B2(0) M"\Bi(0) 

^5 finite. The definition o/C;^'y(R" \ R*') is similar; we need only replace Bj{0) by the tube of 
radius j around R'^, j = 1,2. 
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Although these Holder spaces are our primary tools, we shall also need to refer on occasion 
to a family of weighted spaces and their associated Sobolev spaces. These will be needed only 
for nonisolated singularities, so we restrict ourselves to that case. 

Definition 4 The weighted space r^L'^ on M" \ M!' is defined by 

r^L^{W\m'';r'^-^drd0dy) ^ {w E L^^^(M" \M'=) : J \w\^ r^-^-^^ dr dO dy < oo}. 

The space p^L'^{Q\'E) (relative to standard Euclidean volume measure) is defined similarly, using 
the function p above. 

There are associated weighted Sobolev spaces r^H^{W"' \ M*^) when s is a positive integer, where 
the subscript e signifies that these spaces are defined with respect to differentiations by the vector 
fields rdr, rdy and de. (When s is an arbitrary real number, they may be defined by duality and 
interpolation.) Note that 

C'^'^{W'\R'') Cr^L^{W'\R'') provided cr + (TV - 2)/2 < (5 < 7 + (iV - 2)/2 

C';'°'{n\j:) C p^L^{n\Y.) provided S <-f+{N -2)/2. (11) 



4 Construction of approximate solutions 

Now that the family of radial solutions Ue{x) to (||) on \ {0} has been introduced, we can 
construct the approximate solutions for these problems. We also derive estimates of how far these 
solutions differ from exact solutions in terms of the weighted Holder norms of the last section. 



4.1 Approximate solutions with isolated singularities 

Approximate solutions for (|l|) which are singular precisely at the points of S = {xi, . . . ,xk} 
are constructed by superimposing appropriately translated and dilated copies of the singular 
radial solutions. The K free parameters in the exact solution correspond to the K different 
dilation parameters we can prescribe independently at each of the singular points. First, choose 
X{x) e C^(M^) with X = 1 when \x\ < 1 and x = when |a:| > 2, set xb.{x) = x{x/R). Also 
choose i? > with R < Rq = mii^j (dist{xi,Xj))/2). Let e — (ei, ■ • ■ , } be a K-tuple of dilation 
parameters. Now define 

K K 2_ x-X 

Ue{x) = ^xr{x ~ Xi)u^^{x - Xi) = ^x( — ""'uii -). (12) 

Set /g = Am J + uf . After some computation, we get 

K K 
A = ^ [X - Xi)AxR{x - Xi) + 2\/Ue^{x - Xi)Vxr{x - Xi) + ^(Xr - Xb){x - Xi)uP. {X - Xi). 
i=l i=l 

Lemma 3 For any 7 G M, there exists a constant c, depending on R and 7, such that 
|/c||o,a,7-2 < cEq ^' provided each ei < Cq < 1. 



Proof : For \x ~ Xi\ G (i?, 2R), we get 

\u,{x)\<ce, "-^ 

and 

|Vu,(x)|<ceo 

for some constant c > depending on N and R. The result follows at once. □ 
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4.2 Approximate solutions in the general case 

Now suppose S = ufLi'Si where each is a smooth submanifold in fl of dimension fc^. To 
simphfy the notation here, we assume that aU ki > 0. 

In terms of the Fermi coordinates introduced in §3, the Euchdean Laplacian on T^^*"* can be 
written locaUy in terms of the Laplacians for and Ny'Sf. 

A = Ajv + As, +ei • V2 + e2- V, (13) 

where ei and 62 satisfy 

||ei||o,Q,i + ||e2||o,a,o < c, 

for some constant c > which does not depend on a, nor on x, y. This is discussed in ||^ and 

Now choose a smooth cut-off function x on as before, which only depends on |a::|, and such 
that xix) = for > 2 and xi^) = 1 for < 1. Also, set XRi^) = xi^/R)- If < < 1 and 
R < a/2, then define, in some neighborhood of yo G S^, the function 

2 

Since this function only depends on the variable in the normal space, and is independent of the 
angular variable 9, it is clear that it may be defined globally on all of T^^^ . Now let Ug = J2^i ■ 
As before, let 

Using (p^, we compute that 

fi = ei{x, y) ■ V^Me + 62(2;, y) ■ Wug + u^AnXr + 2Vxi^Vuj + [x^^ - xr)u^k- 



Lemma 4 There exists some c > depending on 7 hut independent of eo < 1 such that if each 
Ci < Co, then 

||/c||o,a,7-2 < CEq, 

where q = — 7. 

Proof : The estimate follows at once from similar estimates already used in Lemma ^ and also 
from the estimates given in ( [l3| ) . □ 
This exponent q is strictly positive by our assumptions on p, provided 7 > . 



5 Outline of the proof 

In both the cases of solutions with isolated or more general singularities, we wish to solve the 
equation (|l|) by perturbing the approximate solutions to this problem we obtained in the 
preceding section. That is, we wish to represent the exact solution u as a sum u — Ug + v, where 
V is small compared to Ug. Thus we wish to solve 

A{ug + v) + {ug + v)P = 0, 

or, what is the same thing, 

Lgv + fg + Qiv) ^ 0. (14) 

Here 

Q{v) = {ug + v)P - - puPg-\ (15) 

is the remainder term, which is quadratically small if v is smaller than Ug in an appropriate Holder 
norm, and 

Lg^A+pur^ (16) 
is the linearization about this approximate solution. 
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We shall show that there is some Cq > and some constant < c < 1 such that for any 
£ ^ coj if ce < < e for all i, there is an exact solution of (|l|) which is a small perturbation of 
Ue- Writing u — Ug + w as before, where w € C^'"(ri \ S), we can ensure that u is singular at 
S provided i' > —"^/ip — 1). In addition, we shall be able to prove that the solution space is a 
if-dimensional manifold. 

When any of the fci > 0, then we shall show later that the solution space is infinite dimensional. 

The analysis of this linearization is the fundamental issue of this paper. It is easy to see that 

L,:C^^"(17\S]) ^e"2(r!\S) (17) 



is bounded for any 7 G M and e. We must show that ( |17| ) is surjective for some 7 > — 2/(p — 1), 
provided that each is sufficiently small. Once this is shown, then it follows that Lg has a 
bounded right inverse 

G,-:C°L"2(f7\S)— .C2."(f7\S). (18) 

We must also show that the norm of this map is bounded independently of the for < eg- 
However, for the range of values of ly for which (^7|) is surjective, it is not injective, so that Gg 
is not unique. As usual, a good choice for a right inverse is the one which maps into a fixed 
complement of the nuUspace of Lg, for example, the orthogonal complement with respect to some 
weighted structure. 

To prove surjectivity of (^7|) we prove both that its range is closed and that the cokernel is 
trivial. When S is a discrete set, the closedness of the range is rather elementary, but when some 
component of E has positive dimension this fact is somewhat deeper. Fortunately, the construction 
of pseudodifferential right parametrices for the general class of elliptic 'edge operators,' of which 
Lg is a particular example, falls within the scope of the theory developed by the first author 
in Existence of such a right parametrix with compact remainder implies that the range of 
Lg is closed in all cases, and that its cokernel is at most finite dimensional for a certain range 
of values of ly. To show that this cokernel is trivial, and that the norm of Gg is uniformly 
bounded for sufficiently small e^, we employ scaling arguments. These arguments proceed by 
contradiction, showing that if first surjectivity, and secondly uniform surjectivity were to fail, 
then counterexamples for increasingly small could be rescaled to obtain some element of the 
nullspace of the global operator Li on M.^ \ {0}, or Li on R" \ M'^ which we show cannot exist. 
Once these results are established, the rest of the proof of the existence of solutions follows from 
a standard fixed point argument. 

We remark here that we have proved these results concerning the uniform surjectivity of Lg 
in slightly greater generality than is used in the nonlinear analysis. In particular, we are able to 
show that Lg is uniformly surjective provided all < eo for eg sufficiently small. However, we are 
only able to prove the fixed point theorem if ae < < e for some e < coi smd for some a G (0, 1), 
i.e. when all the are mutually relatively bounded. 

6 The globalized linearization 

In this section we analyze the behaviour of the operator 

on \ {0}, and later the corresponding induced operator on M" \ M*^. In polar coordinates, Li 
takes the form 

92 TV - 1 a 1 , VJr) 

where Vp{r) — ■ pui{r)P^^ . From Proposition |l| above we have 

lim_o V,ir) = pvP^^ = JiL (at _ ^) ^ A„ (20) 
yj,(r) - cr(2-^)(p-i)+2 as r -> 00. (21) 

Notice that the exponent (2 — N){p — 1) + 2 is negative precisely when p > N/ {N — 2). 
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6.1 Indicial roots 



Li has a regular singularity at r = and, because Vp{r) tends to as r tends to infinity, it also 
has one at r = oo. Hence the asymptotic behaviour of solutions to Liw — are determined by 
the indicial roots of this operator at these points. At r = 0, the indicial roots for Li are 



, UN -2^ ^ 



Here the Xj are the eigenvalues of the Laplacian on 5"^^^, counted with their multiplicity. More 
precisely 

Ao = 0, = iV - 1, j = 1, . . . N, etc. 

The indicial roots for Li at r = oo are the same as for the Laplacian itself, since V tends to at 
oo; these values are 



7f = ^^±y(^^^j J = 0,1,2,.... (23) 

Notice that these numbers are integral and assume all integer values except —1, —2, . . . ,3 ~ N. 
Because Li has a regular singularity at and oo, its mapping properties are well-known, cf. 

Proposition 2 The bounded linear map 

ii : C^;^,(K^ \ {0}) ^ c!;L"2.y_2(R^ \ {0}) 

is Fredholm provided 7 ^ {7^*^ } and 7' ^ {'^^ } . 

We will let 7' — —1 or since we are only interested in solutions which decay at 00, or at 
least are bounded. We could replace this —1 by any value of 7' in (2 — N,Q). However, when 
= 3, we would need to choose 7' G (—1,0). We shall not comment on this further, and persist 
in letting 7' = —1, with this understanding. 

6.2 Numerology 

We record some facts about these indicial roots and some of the other constants that arise 
frequently for later reference. 

Lemma 5 Let p be any number in the range N/ {N — 2) < p < {N + 2)/{N — 2). 

1. The functions —2/ {p — 1) and Ap are monotone in p and 

-2 2-N , 2p f 2p \ N^ - A 

2~N> 7 > and < = — ^ iV < . 

p - 1 2 ^ p - 1 \ p-lj 4 

2. There is a number p* in (j^zg' such that for N/{N — 2) < p < p* , the indicial roots 
7^ are real, with -2/{p - I) < % < {2 - N) /2 < 7+, while if p* < p < {N + 2)/{N - 2), 
then 71^ are both complex with real part (2 — N)/2. 



3. 



"fi < — ^jzTf, in fact ^-^ = ^ 1- Hence in particular, 7j^ are always real. 
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6.3 Injectivity of Li on Cjo(K^ \ {0}) 
Fix two weights i/ and fj. with 

^ < < 5R(7o-) < < SR(7o+) < M (24) 

and jjL + V = 2 — N . These weights are 'dual' in a sense to be explained later. The interval in 
which /i lies is chosen to guarantee injectivity of Li on C^^q; v is determined by this duality, but 
also chosen not to be too negative in accordance with the nonlinear problem. 

Proposition 3 The only solution w e C^'q(]R^ \ {0}) of the equation L\w = Q is the trivial 
solution w = Q. 



Proof: Let {</'i(^),Aj} be the eigenfunctions and eigenvalues for Ag, the Laplacian on ^. 
Then if w is any solution of Liw = 0, it decomposes into the infinite sum 

w{r,0) =J2^,{r)<j)j{0), 

3=1 

where the Wj solve 

N -1 d 1 
= Q^'^i + —^—Q^'^i + ;:2 (^pW - = 0. 

The value j = is omitted in the sum above since every solution of Li,oWo = blows up faster 
than near r = 0. 

Now, by the growth restrictions on w near and infinity, each Wj must be bounded by Cjr'^i 
as r ^ 0, and must decay like r'^~^~i as r — > oo. We shall show that this forces Wj to vanish. 

First of all, notice that if Xj > supp^ then solutions of LijWj — satisfy the maximum 
principle. For these values of j it is also true that 7^ > 0, and so Wj tends to both at and 
infinity, hence must vanish identically. Thus there are only finitely many values of j for which wj 
may possibly not vanish. To deal with most of these, wc use an integral estimate which will show 
that Wj = for Xj > N — 1. Multiply the equation LijWj = by r^~^w and integrate from k 
to R. Integrate by parts and use that > {2 — N)/2 and that Wj decays at least as fast as r^~^ 
as r ^ 00 to sec that the boundary terms converge to zero as /t — > and R — > 00, and that the 
integral converges. The result is that 

/■OO 

/ -{drWjfr''-^ + {Vp{r) - Xj)w']r'^-^ dr = 0. 
io 

Now rearrange terms and estimate V by its supremum: 

/•OO POO 

/ r^-V.«;,f rfr = / {V^{r) - X 
Jo Jo 

/>oo 

< (sup Vp -Xj) r^-^w^j dr. 
Jo 

To proceed further, we note that there is a general inequality for any function which decays like 
r'' at and like r^~^ at infinity: 

rr''-'\wfdr<—^ rr^'-'ldrwfdr. 
Jo - ^) Jo 

To prove this, observe that 

"^~'^wdrW dr, 



,;r^-'\wj\'dr 



pOO 1 poo r\ PC 

Jo N-2J0 N-2J0 



11 



where again the boundary terms vanish in the integration by parts because of the assumed rates 
of decay. Now use the Cauchy Schwarz inequality to prove the claim. 
Use this inequality in the estimate for Wj above to get 



2 roo /•oo 

- / r^-^w^j dr < {supVp - Xj) / r^^^w^j dr. 
Jo Jo 

^2 , 



If supV^ ^ ^ (-^ ^ 2) /4, then this inequality shows that Wj must vanish identically By 
Proposition 0) and Lemma || 

sup Vp < — N 



2 p-1 \ P- 1, 

The function of p on the right is dominated by its value at p = {N + 2)/{N — 2), so that 

N N^-i N(N + 2) 

But now 

4 4 2 

and this occurs for all j except j — 0,1, . . . ,N. 

Thus we have shown that Wj = for j > N, and because of the restrictions on the growth 
at r = 0, we also know that wq — 0. It remains to show that Wj ~ 0, j = 1, . . . N. It turns out 
that we can explicitly write down an explicit solution of Lijw = for this range of j. In fact, 
differentiate the equation Aui + = with respect to d/dxj to get 

Lr|^ = 0. 
dxj 

Since ui depends only on r, dui/dxj — (dui/dr)9j , where 6j — Xj/\x\. Since —AgOj — XjOj, 
j = 1,...N, u[{r) solves Lijw — 0. But u'i{r) decays like r^~^ as r ^ cx) and blows up 
like r^^/'^P^^^^^ as r — > 0. If wj solves LijWj = and decays like r^~^ at infinity and like 

at 0, then some nontrivial linear combination of u'l and wj must decay faster than r^~^ at 
infinity; since the singularities of these functions at cannot cancel, this linear combination is 
nonvanishing. This is a contradiction, since no solution of Cij can decay faster than r^^^ at 
infinity. Hence Wj — 0, and the proof is complete. □ 

6.4 Injectivity of Li on Cjo (K" \ K'^) 

The argument that the induced operator Li = Li + on R" \R'^ is injective on functions which 
blow up no faster than r'^ and decay like r^^ as r — > oo is somewhat more complicated, but rests 
on the results of the last section. 

Proposition 4 The only solution w G Cj^'p (R" \ R*^) which satisfies Liw ^ is w — 0. 

Proof: As in the previous subsection, we analyze this equation by reducing it to a family of ODE's 
and studying these separately. This is done by first taking the eigenfunction decomposition in 6, 
as before, but then also taking the Fourier transform in y. Letting rj be the variable dual to y in 
the Fourier transform, we obtain the family of operators 



|2 



If Liw = and w S C^'q i then 

oo ^ 

w{r,e,y) j e'y'' Wj{r,r])dy, 
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where Lj-j^pw^ = and wj € C^'g (R+) and depends distributionally on rj. Clearly, if L^j^piZij = 
with lij G C^'o (R+) imphes Wj = for every j, then Liw = with w G C^'g (R" \ M'"') implies 
w = as well. 

First observe that the integration by parts argument of the last subsection may be repeated 
essentially verbatim to show that all solutions of this equation are trivial unless j = 0, 1, . . . , iV. 
Furthermore, since /i > 5R(7^), there are no solutions in C^'g (M+) for j = 0, either. Thus we are 
left to consider only the case j = 1, . . . N, where Xj = N — 1. 

For simplicity, in this argument, let us drop the index j in our notation and replace by 
a parameter E. Then Lj |,j|2 becomes L^;, and so forth. Now let us recall that any solution of 
ILew = has a Frobenius series expansion (convergent in some interval, since Vp is real analytic) 
of the form 

oo oo 

w^J2 a^r'^'^' + J2 bjr''^^^- (25) 
j=o j=o 

There are additional terms, with logarithmic factors, in the second series in the case when the 
two indicial roots differ by an integer. Since the argument is quite analogous in that case, we only 
consider the case when these extra terms do not appear. The two leading coefficients, oq and bo 
are free, and may be specified arbitrarily; all the aj for j > are determined in terms of ao by a 
recursion formula, and similarly for the bj. In particular, if ao = 0, then all Uj = 0. Hence, since 
^ > 7]^, we have that w G C^'" near r = if and only if ao = 0. 

We shall show by a sort of continuity argument, treating both p and E as parameters, that if 
w is a nontrivial solution on R+ growing at most polynomially as r — > cxd then ao never vanishes. 

For any two numbers < ro < ri < cx) there is a map Crg.n which sends the Cauchy data 
(M;(ro), ^'(ro)) of a solution w at ro to its Cauchy data {w{ri),'w'{ri)) at ri. When tq = 0, this 
Cauchy data is replaced by the pair of leading coefficients (ao,6o) in (|25|). If < rg then it is 
well-known that Crg^n is invertible and depends smoothly on the parameters {p,E). It is also 
true that if ri is in the interval of convergence of the series ( p5| ) then Co,ri is invertible and still 
depends smoothly on (p, E) . 

Amongst the two dimensional family of solutions of L^w = 0, there is one solution, unique up 
to multiplication, which decays at infinity. Actually, if "hsw — and if w is of polynomial growth 
as r ^ oo, then w decays like r'^e^'"^ for q determined by the coefficients of L^. We fix the 
solution uniquely by requiring that w{\) = 1. This makes sense because, by the normalization in 
the function Vp{r) is less than some fixed constant, in particular less than — 1, for r > 1. 
Thus solutions oI'Lew = satisfy the maximum principle on [l,oo). Hence any solution which 
vanishes at r = 1 cannot decay as r — > oo. 

For this decaying solution w look at the leading asymptotic coefficient ao — og (p, E) in its 
series expansion at r = 0. We write out this dependence on the two parameters p and E explicitly 
since we need to show that ao(p, E) ^ G for (p, E) in the strip (A^/(Af-2), {N+2)/{N-2)) x [0, oo). 
For if this holds, then w ^ C^;o (R+) and our result will follow. 

We first show that ag is smooth in {p,E). We do this by using the maps Crg^n discussed 
above. First consider the Cauchy data (1, w'(l)) of the decaying solution w. Since the coefficient 
of order zero of L^; is negative on [l,oo), w may be constructed by the 'shooting method' there, 
i.e. by solving the boundary problem L^w^ = on with w/j(l) = 1 and Wfi{R) = 0, 

and letting R — ^ oo. By the maximum principle, the limit must exist, and since it is bounded 
it agrees with w. Since the Cauchy data (l,w^(l)) depends smoothly on {p,E) for each R, it is 
not difficult to see that the limiting values (l,u;'(l)) do as well. Now use the composition of the 
(invertible) maps Co,ro for ro very small, and Cr„^i to conclude that (ao,feo) depends smoothly 
on ip,E). 

Now we proceed with the rest of the argument. First observe that ao(p, 0) ^ for any p 
by the result of the previous subsection. In fact, writing out the explicit solution dui/dr more 
explicitly, we can check that ao(p, 0) > 0. Next, note that for p sufficiently close to its lower limit 
N/{N — 2), the function V stays uniformly small, and in particular less than A^ — 1 for all r. 
Hence we can again apply the integration by parts argument of the previous subsection, but now 
using that sup{Vp — A^ + 1) < 0, which immediately shows that ag 7^ for all p < N/ {N — 2) + n 
for every E and for some k > 0. By continuity from £^ = 0, ag > for (p, E) in this smaller strip. 
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Finally, we can apply a similar argument when E is sufRciently large, regardless of the value of p. 
In fact, all we need is that the suprcmum of Vp{r) — N + 1 — r^E is less than {N — 2)^/4, and this 
is clearly true for any p, for E large enough. Thus, ao > also for (p, E) in this range of values. 

Finally, suppose that there is some value of (p, E), with p G {N/{N - 2), (TV + 2)/(iV - 2)), 
for which oq = 0. Using the regularity of ao in these parameters, there is some {p^, Eq) for which 
ao(po,£'o) = and (ao)_B(po7 -E-o) = (i.e. both aq and its derivative with respect to E vanish 
there). Differentiate the equation L,ew = and set {p, E) — {po, Eo). li w — dw/dE then 

Since ao = 0, the right side of this equation blows up only like r^^ as r — > 0. Thus w is the 
sum of two terms; the first blows up no faster than r"' while the second is a solution of the 
homogeneous equation and might blow up like r'^ . However, since the leading coefficient of 
w is the derivative of the leading coefficient of w, i.e. it is (ao)_E, which vanishes, we see that 
this solution of the homogeneous equation is absent and w blows up no faster than r''^"''^. Now 
multiply the equation by w ■ r^~^ and integrate from to infinity. Using the exponential decay 
as r tends to infinity, and that w ^ r'~' , w ^ r'^ +^ as r tends to 0, we see that this integral is 
well-defined, since 27+ > 2 — N. Furthermore, it is permissible to integrate by parts to get 

/•OO /"OO /"OO 

0=/ {LEw)wr^-^dr = w{hEw)r^-^ dr ^ 2 \w\'^r^^^ dr. 
Jo Jo Jo 

Hence w = for this value of (p, E), which is a contradiction. □ 

6.5 Surjectivity of Li and Li on C^'"i 

From the results of the last subsection we may deduce the following: 

Proposition 5 Let ^^j^ <v = 2 — N — 5R(7(7) as before. Then the maps 

L, : Cl:^,iR^ \ m ^ C^"2,_3(R^ \ {0}) 

and 

Li : C^;!;i(R" \ R*^) — > C|!:"2,-3(K" \ 1^'') 

are surjective. Furthermore, the first map is Fredholm and its nullspace is one dimensional, while 
the second map has an infinite dimensional nullspace. 



Proof: We shall only sketch the proof of this result and provide references to papers where 
analogous facts are proved thoroughly. As discussed already in §5, the proof of this result is much 
simpler for the former of the two maps. In fact it is not difficult to write down a right inverse for 
Li directly in terms of solutions of the homogeneous problem Liw = 0. In order to do this, it is 
crucial to use the fact that any solution which is bounded by a multiple of r^^ near r — does 
not decay as r tends to infinity, and any solution which decays at infinity must decay at least as 
fast as r^~^. Note that this is equivalent to the injectivity of Li on C^li. The right inverse for 
Li may be written as a sum of the right inverses for Li ^ on each eigencomponent; each of these 

right inverses may be constructed explicitly from the solution which blows up like r'^'j as r — > 
which is unbounded as ?' — > 00 and the solution blowing up like r^j as r — > and decaying as 
r — > 00. It is then not difficult to show that this sum of right inverses for each eigencomponent is 
bounded on the weighted Holder spaces, cf. for this procedure for a closely related operator. 
This proves that Li on C^'"]^ is not only Fredholm but surjective. 

Any solution of Liw = in C^'li must be radial, by considering its growth as r — > 0. There 
is a two dimensional space of solutions of Li^wo = 0. By considering a suitable nontrivial linear 
combination of any two basis elements of this space we obtain a solution which decays like r^^^ 
as r — > 00; there cannot be two independent solutions with this decay rate for the same reasons 
as in the previous subsection. This solution grows like some combination of r^o as r — > 0, but in 
either case is in C^'li. Hence the nullspace is one dimensional. 
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To establish the corresponding facts for Li we use Hilbert space techniques and the construc- 
tion of a pseudodiffercntial right parametrix for Li from Li is surjective on C^li provided 
its range is both closed and dense. The existence and boundedness of this right parametrix gives 
the closedness of the range. Then duality, coupled with Proposition ^, yields that the cokernel of 
this map is trivial, so that the right parametrix may be replaced by a right inverse. 

It is most natural to construct the pseudodifferential right parametrix for Li relative to the 
spaces r^L^ because of the central role of the Fourier transform in this construction. Once the 
parametrix is obtained, it is then necessary to show its boundedness on the weighted Holder 
spaces. These steps are carried out in detail in 1^, to which we refer the reader, cf. also |8|. Here 
we simply show that Li satisfies the hypotheses necessary for that machinery to apply. 

The main result of |^ implies that Li has closed range on r^L^ provided S ^ i^f}^ where 
5^ = + {N — 2)/2. Choose a number —Si, just slightly smaller than + {N — 2)/2, in 
particular so that -2/(p - 1) + {N - 2)/2 < -S„ < ly + {N ~ 2)/2 < 0. Then -6„ ^ {6f} and 
C^'"(S1\I]) C It is also proved in that Li is essentially surjective on r~^" L^, i.e. 

it has closed range and at most a finite dimensional cokernel, and that it is essentially injective, 
i.e. has closed range and at most a finite dimensional kernel, on r^^"L^. These two facts are 
consequences of one another by duality since r*L^ is the dual space of r~^L^ and Li is self-adjoint 
on r"!/^. The crucial hypothesis that must be satisfied in order for these last two conclusions 
to be true involves the 'model Bessel operator' Li for Li. This operator is obtained by freezing 
coefficients of Li at r = in an appropriate sense, taking the Fourier transform in y, and finally 
rescaling by setting s — r\r]\. The operator obtained in this way is 

L _ ^ d ^ Ap + Ae ^ 

ds^ s ds 

This crucial hypothesis is that Li is injective as a map on s''''L^(M+ x S'^^^; s'^"^ ds d9), or 
equivalently, surjective as a map on s^^'' L'^(R'^ x S^~^; s^~^ dsdO). 

Li may be analyzed directly to show that this hypothesis is satisfied. In fact, introducing the 
eigenfunction expansion with respect to Ag , as usual, we obtain a family of ordinary differential 
operators, the solutions of which may be determined explicitly in terms of Bessel functions. 
Analogously to the situation for Li we find that any solution which blows up no faster than r'' 
as r — > grows exponentially as r — > oo. Using these solutions we can construct a right inverse 
for Li explicitly. 

Since Li satisfies the hypothesis, we conclude that Li itself has closed range, with at most a 
finite dimensional cokernel, as a map on Cl'"^. The cokernel of this map may be identified with 
the kernel of Li as a map on r^" L'^ or C^'"]^. But from the previous subsection, we know that 
this cokernel is trivial; hence Li is surjective, as desired. 

Finally, we may explicitly exhibit an infinite dimensional nuUspace of Li in Cl'li- For, if 
Liw = and w is bounded by as r — > 0, then only the eigencomponent wq may be nonvanishing, 
just as for Li. Now, taking the Fourier transform in y, we see that Lq i^pt&o — 0. Hence zi'o(''j v) = 
A{rf)WQ{r,ri), where Wq is the unique solution of this equation which decays (exponentially if 
7y 7^ 0) as r ^ 0, normalized so that Wo(l, ??) = 1- The coefficient Airf) is allowed to be arbitrary, 
so long as the corresponding is bounded uniformly, for each r, as \y\ — > oo. Clearly there is 
an infinite dimensional freedom in choosing such coefficients; for example, we could let A{ni) be 
an arbitrary element of the Schwartz space. This completes the proof. □ 

7 Injectivity of on C^Jl^ \ S) 

We are now able to turn to the first of our main tasks, to show that is surjective on C^p(f2\I]) 
when all the are sufficiently small. By an argument identical to the one indicated in the proof of 
Proposition 1^, this surjectivity is equivalent to the injectivity of this operator acting on functions 
growing like at the singular points, where (2 — N)/2 < /x = 2 — iV — j/ < 7^ < 0, as in the 
previous section. This injectivity is what we shall prove. 
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7.1 Preparatory Lemmas 

There are three lemmas which are used in the rescahng proofs in both cases. We shah adopt the 
notation 

Here, using shghtly different notation than in previous sections, B^.{T,i) is the tubular neighbour- 
hood around E.^ of radius e^, and, in case S is discrete, E.^ simply equals the point Xi. 
The first Lemma states that Lg satisfies the maximum principle on fig. 

Lemma 6 There exists some eo > such that if all Ci < eo, then after a normalization of the 
initial radial solution ui, the operator Lg satisfies the maximum principle on D,g. That is, if w < 
on dflg and LgW > 0, then w < on Qg. 

Proof: We give the proof of this result in the case where E is discrete. The general case, where 
E has positive dimension can be treated similarly. Set = maxjit;, 0} on fig and w'^ = on 
each (Ei), so that w'^ G i/o(17). Multiply the inequality Lgw > by w'^ and integrate by 
parts to get 

K 



I \Vw+\^< f \w+\^puPg-^ <paP-' f (V 



x-xA-^)\w+\^. 



Here a is the supremum of wi = r p-i ui for r > 1, which by Remark y can be taken as small as 
desired. Now, using the identity 

(iV-2)|a:|-2 = div (^), 
and integrating by parts, we obtain the inequality 

If a is taken sufficiently small so that the inequality Apa^^^K < {N — 2)^ holds, we see that 
w'^ = 0, and hence w < 0, as desired. □ 

The second Lemma uses this maximum principle to deduce estimates in the weighted Holder 
spaces for solutions of LgWg = fg in terms of their boundary values on dflg. 

Lemma 7 There exists eo > such that, if all ci < sq then the following estimate holds. If 
fg G C^^2(^e); where 7 G (2 — A'^, 0) is fixed, and suppose that Wg is any solution to LgWg = fg 
with Wg — on the 'outer boundary' of fig, i.e. on dflg Pi dfl. Then there exists a constant c > 
independent of e such that 

\\wg\\2,a,y < C ^||A||0,a,7-2 +^e7'^lhf|lo,0,aB,.(Ei)j , 

where the first two norms are taken over Qg. 



Proof : Define 4> £ C^{fl \ E) to be a positive smooth function for which, in some fixed 
neighbourhood i?o-(Ei) for each i, (j){x) = dist {x, E^)"*"; here a should also be chosen so that the 
approximate solution Ug is supported in Ui3o-(Ei). If some E, is a point, Ej = Xi, then — [x — Xil"' 
in Ba-{xi). For example, we can take </> = p'*', where p was the function introduced in §3. Then a 
simple computation shows that in B„{xi) \ B^.(xi) we have 

L,Jx-Xip = \^-f{N -2 + J)\x-x^p-^ +puPg'^\x-x^p^ < -c\x~x,y~^, 

where the constant c > can be chosen independent of e, since 2 — < 7 < implies that 
7(iV — 2 + 7) < and since pa^"-'^ can be chosen as small as desired. If E^ is positive dimensional, 
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then a similar estimate holds. For this we simply need to use the expression (|l^) and the estimates 
for ei and 62 there. 

Let A > denote the supremum of on lJ^idBc{xi). Hence (p, multiplied by a suitable 
constant times (^A + | |/e| |o,a,7-2 + X^ili ^^'''ll^cllo.o.as,. (so) ^ supersolution for the problem 
in B„{xi) \ B^.{xi). Likewise, —0 multiplied by a similar constant is a subsolution. 

We claim that, if eg > is small enough, then A is bounded by a constant times 



^||/£||0,a,7-2 +^ei ''l|u'e||o,0,aS,,(S,)^ ■ 

In order to prove this claim we will argue by contradiction and assume that we have a sequence 
of counterexamples to this assertion. Thus suppose that there is some sequence of i^T-tuples 
e(^) = (ef ^ . . . , e^-*) such that for some subset of the indices I, . . . , K, the corresponding e^p tend 

to zero. For convenience, we take this subset to be {1, . . . , J}, so that e''p — ^ for j < J and 

(£) . . (£) 

> c > for j > J. In addition, up to a subsequence we may assume that converges to 

some ej > c for all j > J. Suppose also that for each such e^^^ there is a function w^a) G C^'"(rig(*) ) 

with LgWgit) = fg(i) and such that w^(e} = on the 'outer boundary' of fig and fg(i) G Cl^^2(^e) 
bounded independently of £. First multiply the w^ii) by a suitable constant so that 

sup \wg(t) I — 1. (26) 

Using the super-solution constructed above, it is easy to see that Wg(t) converges to a solution w 
of Aw = in \ u/^j^Ei 'jfLj+i B^ii^i)- In addition, w also vanishes on dfl ufLj^i dB^-CSi), 
and finally w e C^'". Since 7 > 2 — iV it is well known that the singularities of w at u/^j^Ei are 
removable. Therefore w solves Aw = in \ U^j^j^i?e; (S^) with boundary data, so w = Q 
which contradicts (^6|). The claim follows. 

We conclude that there is a constant 6 > independent of e such that Li{wg + 5(f) < fg — cS(t) < 
and Lg{wi — 5(f) > fg + cScj) > 0. By increasing 6 to make Wg — 6(j) < and Wg + 5(j) > on dilg, 
we deduce that \wg\ < 6(j) in all of fig, i.e. w^ e Cj^'°(rie), with norm independent of e. Standard 
arguments using rescaled elliptic estimates show that Wg G C'^'°'{flg) with norm independent of e. 
□ 

The third Lemma shows that the weighted Holder norm with exponent fi for a solution of the 
fixed homogeneous equation Liw = on a fixed tubular neighbourhood B^CSi) is controlled by 
the norm of w on dB^CSi). This is quite obvious when E^ is a point, since the nuUspace of Li is 
finite dimensional then. When E.^ has positive dimension, this is a more substantial result, and 
it follows directly from the results of [^. Nonetheless, we give an elementary proof, which covers 
both cases. For convenience here we drop the subscript i. As usual, (r, 9, y) are Fermi coordinates 
on this tubular neighbourhood. 

Lemma 8 Suppose W£ is a sequence of solutions of LiWg = in B„{T,), with Wg G C^'"(i?cr(E)) 
and \wg\ < A on dB^{^) uniformly in I. Then \\Wi\\2.a,iJL cannot diverge as £ ^ 00. 

Proof: First of all, observe that by the rescaled Schauder estimates, it suffices to show that 
the supremum of r~'^|w£| over B^ does not diverge. Suppose, to the contrary, that it does; 
furthermore, suppose that this supremum takes the value Ci, which tends to infinity, and is 
attained at some point xi G B^. Let (re,6i,yi) be the Fermi coordinates of xg. Since wg can 
be bounded by some new constant A' on i?a(E) \ B„^2{'^)i the rg must converge to zero. For 
if they did not, then consider the rescaled function wg = C^^wg. This still solves Liwg = 0. 
Furthermore, the supremum of r~^'Wg on B^[Y?) is equal to one, and is attained in some fixed 
annulus B„{T,) \ B[}{T,). Thus, the wg converge to a nonzero limit Woci which is a solution of 
LiWco = in _Bct(E). However, the supremum of wg on B^ \ B^/2 tends to zero, so that iDoo = 
on this annulus, which is a contradiction. 

Now, since we have established that rg — > 0, we may again rescale. If E = is a point 
(which is the origin in these coordinates), consider the function wg{r,9) = C'[^r^wg{r /rg,6). If 
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E is higher dimensional, then some subsequence of the points yi converge to ?/oo G ^ (since S 
is compact). Choose Fermi coordinates centred around this point, so that j/oo — 0, and these 
coordinates are defined for \y\ < t. In this case let wg — C^^r'^'We{r/rg,9, {y — y^)/ri). In the 
former case, wt satisfies Li/j-^we = on i3g./j.j(0), and is bounded by there. The same is true 
in the latter case, except that the operator Li/^^ must be replaced by one for which the error 
terms ei, 62 are replaced by very small translates, and which are in any case still tending to zero 
when £ ^ 00. This rescaling has been chosen so that r^'^\wi\ attains its supremum on dBi{T,). 
As before, pass to a limit, Woo- By the previous remark, Woo 7^ 0. 

When E = 0, then it is also true that Woo G Cf^''^{lS.^ \ {0}), and is a solution of 

there. But the solutions of this equation are of the form 

00 

j=o 

It is clear that no function of this form can be bounded by both as r — > and r —>■ 00. Hence 
we arrive at a contradiction again, so the assertion of the lemma must be true when E is a point. 

In case dim E = fc > 0, we may still take a limit, and get a function Woc which solves the same 
equation, but now on R" \ ]S.^. It also satisfies the estimate supr~'^|?Z)oo| < 00. However, intro- 
ducing a decomposition into eigenfunctions for the spherical Laplacian Ag yields the uncoupled 
system of equations 

f 9^ N~l , , „ 1 , X 1 A 

[q^ + —;— + + ip< - ^1) ^) ^1 = 0' 

where we have renamed the eigencomponents simply Wj. Each such eigencomponent satisfies 
supr^^lwjl < 00 uniformly as \y\ — > 00. Taking the Fourier transform of this equation in y 
reduces it to an equation of Bessel type: 

This equation may be solved explicitly, and as discussed earlier it is easy to see that for 7^ the 
only solution of this equation which grows no faster than as r — > must increase exponentially 
as r ^ 00. Thus Wj{r, 77) = for r] 0, and so it is polynomial, and hence constant in y (since 
Wj is bounded in y). Therefore, the previous argument applies to show Wj — for all j. Hence 
Woo — 0, which is again a contradiction. This proves the Lemma in all cases. □ 

7.2 Injectivity 

Theorem 5 If eo is sufficiently small, and if each e,; < eo, then 

L,:C^;?,(f^\E)^C;:L"2(f^\S) 

is injective. 

Proof: As in the previous lemma, the proof uses rescaling to argue by contradiction. The 
idea is quite simple. If the result of this theorem were false, there would exist a sequence of 
counterexamples, corresponding to some sequence of the e^^\ with at least some subset of the 

decreasing to zero. By rescaling these functions and passing to a limit, a counterexample to 
Propositions || and ^ would be obtained. 

Thus suppose that there is some sequence of iiT-tuples e^^^ = (e^^\ . . . e^^) such that (possibly 
passing to a subsequence) for some fixed subset of the indices 1,.. . ,K, the corresponding 
tend to zero. For convenience as before, we take this subset to be {1, . . . , J}, so that ej^' 
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for j < J and > c > for j > J. Suppose also that for each such e^^^ there is a function 
We G C^'p(ri \ E) with L^(i)We = 0. First multiply the wt by a suitable constant so that 

sup p{x)-''\wi\ = 1. (27) 

From Lemma 0, we get that 

sup p{x)~^\wi\ < C sup p~^\we\ < C. 
uas^(f)(Si) 

For some subsequence (which we assume is the full sequence), the suprcmum of p(x)^'^\we\ on 
UdB («)(Ei) is attained on some fixed dB (e){'Ej). 

i .... . ^ 

There are several possibilities, which we examine in turn. For simplicity, we initially consider 
the case where S is discrete. First suppose that for the sphere where supremum above is attained, 
the index j satisfies j < J. Fix this index j, and for convenience, translate fl so that Xj = 0. 
Then ej^^ 0. Now rescale, setting 

iS,(a;) = (efV^«Hef2;). 

This function is defined on i? . («)(0) and solves (A + pu^^^)wg = there. Here ui is the radial 

solution ui{r) multiplied by a cutoff function xi^\^^x/ R), which is equal to one on an increasingly 

large set as e^^^ tends to zero. By construction, the suprcmum of r^'^jiD^I is bounded by C on 
B ,,(.,(0)\Bi(0). 

^ n 

We can apply Lemma H to conclude that the r~^|iSf | are bounded uniformly on Bi{0). Now 
let £ — > oo. Since their norms are bounded, the we tend to a limit Woo G C^'"(i?i(0)), which is 
defined with r^^jwool bounded on all of \ {0}. In addition, it is a solution of LiWoo = 
there. But Woo is not identically zero, since its supremum on dBi(0) is one. This contradicts 
Proposition |[ 

This last argument still leads to a contradiction if it is only true that for some j < J (and 
some subsequence of the £), the supremum of p^^\we\ on dB (e)(xj) is bounded below by some 

positive number C > 0, independently of £. 

Having excluded these cases, we can now assume that 

sup p{xy\we\ — > 

every j < J. 

The norms |kL'£||2,Q,^ (over all of f2) are bounded. This is clear in the region flg(e) by (p7|), 
and by Lemma ^ it is also true for the balls with index j > J. For the balls with index j < J, 
note that the rescaled functions iue — {e"p)'''^we{e^p x + Xj) on -Bi(O) have the same weighted 
norm there as the unrcscalcd functions do over B {i) (xj). Since the supremum of iJue over dBi{0) 

is bounded by (e^)^^, which tends to zero as — > oo (and in particular, is bounded), we can 
apply Lemma ^ again to conclude boundedness of the weighted norms of we in these balls too. 

Now (pass to an appropriate subsequence and) let £ tend to infinity. The we converge to a 
nonzero limit Woo, since they have supremum equal to one on dB (f){xj) for some j > J, i.e. on 

a sphere which has radius bounded away from zero. This limit is an element of C^'^(f2 \ S), and 
solves LgiWoo = 0, where e' has e'j = for j < J. Hence the potential pu^~^ in Lgi is singular 
only at the points Xj for j > J. Since Woo G Cf^'" locally near each xj, and since p, > 2 — N, 
it easy to check that Woo is a weak solution of Lgiw = in a full neighbourhood of the points 
Xj with j < J. Since the operator is smooth at these points, a standard removable singularities 
theorem shows that Woo is smooth except at the Xj , j > J. This means that Woo is a solution of 
this operator with singularities at some discrete set E' with strictly fewer elements than E. Now 
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we may proceed by induction, the case K — I already having been treated by the proof above. 
This completes the proof when S is discrete. 

The modifications necessary to handle the general case are rather minor. In fact, it is only 
necessary to modify the way in which the rescaling is done and invoking Proposition ^ at the 
appropriate place. Thus, starting the proof in the same way, if the functions r^^jw^j attain 
their maximum at a point zi — [ri, 9i, yi) in Fermi coordinates around some E^, then we use the 
rescalings 

Mr, 0, y) = iefrwdr/ef,d, (y - yd/tf). 
The rest of the proof is analogous to the previous case. □ 

Remark 4 It is actually possible to prove that in the case of nonisolated singularities, the rescaled 
limiting function ?Doo is independent ofy. This may be done by using some regularity results from 
to show that the tangential oscillation of the initial sequence wg may be bounded independently 
ofe. 



8 Uniform surjectivity of Li on \ S) 

The second main step in the linear analysis is to show that there is some choice of right inverse 
Gg of Lg on C|^^2(^ \ ^) which has norm bounded independently of e provided each ej < eq. As 
we have indicated before, the subtlety here is that Gi is not unique, since is not injective on 

'^y ■ 

The usual choice for Gg is as the right inverse whose range equals the orthogonal complement of 
the nuUspace of L^. Of course, since we are working in Holder spaces, this orthogonal complement 
is meaningless. However, what amounts to the same thing is to require that the range of G^ lies 
in the range of a fixed adjoint of Lg. Again, this adjoint depends on some Hilbert space 
structure, but once we have chosen an appropriate one with respect to which the adjoint is taken, 
we can forget about it and simply use this adjoint. 

By ( pd] ) and using the notation of §6.5, C^'" is contained in p"^" L^{p^~^dr d9 dy). In the 
following, the space will always be taken relative to the background Euclidean measure. We 
shall consider the spaces p^L^ and p^^L^ to be dual with respect to the natural pairing 

p'L^xp-'L^ 3 {^,4,) J H. 

Relative to this pairing, the adjoint of Li : p^^" LP' p~^'^^^L'^ is just : p^-'+^L^ p^"L'^. 
We have proved that the former of these maps is surjective and the latter one is injective. Using 
the fixed isomorphisms 

2(5 ~5t2 St2 

p : p L — > p L , 

we may identify this adjoint, L^, with p^'^^" L^p^^"^'^ . Premultiplying by we can regard 

It is not hard to see that Lemma ^ is valid for both and Lg . 
Now form the fourth order operator 

This map is an isomorphism. Hence there exists a bounded two-sided inverse 

g,:p-'^-^L^^p-'"+^L\ 

In particular, CgQg = I, which means that Gg — L*^Qk is a right inverse to which maps into 
the range of as desired. Henceforth, Gj will always denote this particular right inverse. 
It can be shown directly from the structure of Qi and Gj (see Q) that 

g,:Cl-^^{n\Y.)~^Ctt2,v{^\^) 
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and 



G,:C^,{n\j:)^cl:U^\j:) 



are bounded. In the former of these the subscript P denotes a more general set of eUiptic boundary 
conditions for this fourth order operator, i.e. the domain is restricted to functions u which vanish 
along with L'^u at dfl. 

We can now turn to the proof that the norm of Gg does not blow up as the components of e 
tend to zero. This uniform surjectivity is an immediate consequence of the following two results. 

Proposition 6 Let e be some K-tuple with all e.j < eq. Then provided eq sufficiently small, 
there is no solution of the system of equations LgU = 0, u = L^v with u £ C^p(r2 \ S) and 

e Cl'^2.vi^\^)- 

Proof: This is really just a corollary of Theorem |[ In fact, suppose that u and v satisfy this 
system. Then L^L^v — 0. Recalling that is identified with p~^^'' Lp^^" here, multiply this 
equation by p^^'^v = w and integrate with respect to standard Euclidean measure. Then the 
integrations by parts in 



are valid because d^w G C^^i^'g' ^'^^ ^ — 4' ^^'-^ because w satisfies Dirichlet conditions on 
dft. (Here r is the Fermi polar distance coordinate near each S^. Also note that we can choose p 
to be constant near dil so that p'^^"v and L*^{p'^^"v) still vanish at the boundary.) Hence L^w = 0. 
But since w G C^j" for some p' > K(7^ ) , we may conclude from Theorem ^ that w = provided 
eo is small enough. □ 

Theorem 6 Let e*-^-* be any sequence of K-tuples, with each e'^ < eg. Let f^iQ be any sequence of 
functions in C°^2(^\^) '^^^^ norm uniformly bounded as £ tends to infinity. Let u^it) € p(r2\I]) 
be the unique solution of LgU^(e-) — f^{t) which also satisfies u^a-) € ran {L*^(y), i.e. u^[i) — L*^(yVgie) 
for some Vg{e) € C^!j^2 v Then the norm of u^n) in C^'" is bounded uniformly as £ oo. 



Proof: This argument is again by contradiction, and is very similar to the ones in the last section 

4 



Clearly we only need to consider the case where some subset of the tend to zero. As before, 



we assume that ej^-* for j = 1, . . . , J and ej^^ > c > for j = J + 1, . . . , To start, by 
hypothesis 



supp-'^+^l/.wl <C, 
n 



for all £. Applying Proposition we get 



Now define 



sup p "lugd)] < C + C sup p '^\ug(i)\ 

n\uB (f) (xi) UdB (f) (xi) 



Ag = sup p ''\ug(t)\ 

UdB 



We resize the functions , Vg(£) , fg{e) by setting 

Ugl() = Ug(t) /Ag, = f^fi) = fg(i)/Ag, 

so that Lg(e)Ug(fj — fg(e) and — u^it). If A^ stays bounded as £ tends to infinity, then we 

are finished. If not, as we now assume, then | | |o,a,i/-2 tends to zero. 

We wish to take a limit of the equations Lg(f)Ug(e) = fg{e) and L*^i-^Vg(() = u^it). The main 
point is to show that and tend to limits. For each £, choose a point zg where p~'^Ug{t) 
attains its maximum. By passing to a subsequence, we can assume that zi stays in some fixed 
B it) (Ei). 
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The simpler case is when i > J, so that does not tend to zero. Since by the above and Lemma 
^the supremum of p~'^\u^{t) \ in B_{t) (Ej) cannot diverge and stays bounded away from zero. Then 

we can directly pass to a limit as £ tends to infinity and obtain a function u G C^'^{VL \ S) which 
is nonvanishing and satisfies L^m = 0. Here R is the [K — J)-tuple comprised of the limiting 
values of the for j > J. Since u blows up no faster than p'^ at the Sj with j < J, and 
since the limiting term of order zero in the operator is smooth at these submanifolds, the same 
removable singularities theorem as we used before shows that u must be smooth at those points. 
Next, since the norm of is uniformly bounded, we can also pass to a (weak) limit and obtain 
some function v G Cj^^2i)(^ \ ^) with L\v = u. Since u is nonvanishing, v must also be. In 
addition, by elliptic regularity, v is smooth at the Ej, j < J. However, this is a contradiction to 
Proposition ^ Thus we may assume that for the index i for which the maximum of p~^\Uf(i) \ is 
attained in _B_(f)(Si), the radius e'p tends to zero. 

Let us first treat the case where is a point. Translate the whole problem so that = {0}. 
Next, rescale the two equations by the factor , in particular, u-^n) is replaced by 

and similarly for ^^(f) and j^it). For convenience, replace ef^ by e. Then in the ball i3o-/e(0) 
we have tending to zero in C^^jj '^'^^ '^cW having norm in C^^" bounded uniformly by one. 
Replace the rescaling of the function p now by the polar variable r in this expanding sequence of 
balls. In the new rescaled coordinates, let the supremum of r~^\u^{t) \ be attained at some point, 
which we still call xg,. Then \x(\ is bounded. By the argument used in the previous section, we 
see that \xi\ must also stay bounded away from zero, otherwise we could rescale by the factor 
jx^l"^ and arrive at a contradiction, as in that section. Thus, there is some uniform lower bound 
< C < \xi\ < 1, and hence we can pass to a limit of the Ugii) in C^'". Let us call the limiting 
function u. Then 

Liu = in \ {0}. 

We wish to show that Vg(e) also tends to a limit. To see this, first note, that by the maximum 
principle again, 

sup r~''\vg(t)\ < C + C sup \vg(e)\. 

S„/.(0)\Bi(0) aBi(O) 

Define 

A'^ = sup r"'^|wj(f)|. 

Si(0) 

If A'^ stays bounded as e — > 0, we can take a limit of the v^(e) to get a function v G C^iJ^jl^^^ \ {0}) 
such that Llv — u there. 

So, suppose not, i.e. suppose A'^ tends to infinity. We resize all the functions once again, by 
letting = u^{i)/A'^, and so on. Then = /^(^) and L\v'-^i.^ = u'--^^). By construction, 

the norm of f'-^^^ tends to zero in C|^!^2 ^^"^ the norm of u'^^,^ tends to zero in C^'", but the norm 
of v'-^i) in C^!^2 stays bounded. If the supremum of r^^^'^\v'-(^i) \ occurs at some point xi, then 
exactly the same arguments as above show that \xi\ stays bounded away from and oo. Hence 
once again we could pass to a limit to get a solution V G C^^J^jC^^ \ {0}) such that L\V = 0. 
But we know that L\ has no elements in its nuUspace which decay at infinity, which implies that 
F = 0, a contradiction. 

Finally, then, we have arrived at the situation above, that there exist u G C^'"(M^ \ {0}) 
and V G C;:+2(l^^ \ {0}) such that Lm = 0, LI V = u. Combining these two equations, we get 
LiL\v = 0. But since v decays at infinity, we can multiply this equation by v and integrate. 
Integration by parts now shows that L\v — 0, which we have already observed implies that w = 0, 
hence u = as well. This is the final contradiction. The only alternative is that the C^'" norm 
of can not blow up. 

The case when some are of positive dimension is treated very similarly. In fact, only the 
rescaling needs to be done slightly differently, but in the same manner as at the end of the proof 
of Theorem |^. □ 
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9 The fixed point argument 



We are now in a position to complete most of the proofs of Theorem ^ and Theorem ||; what 
wiU remain after this section is the assertions about the moduh of solutions. We need to find a 
solution G C^'^(ri \ S) to the equation 

L,v + Q{v) + /e = (28) 

where /j is the error term introduced in §4, either for the case when E is discrete, or in the more 
general case, and Q is the quadratic remainder (ffq). We do this by the standard contraction 



mapping argument. We will define a continuous operator K. from the space C1'^{VL \ E) into itself 
and then prove that this map is a contraction on some small ball in this space. 
li V \s a. solution of (ESf) , then + w is a weak solution of 



-A('Ug + i;) = + in Vl 

(29) 

Uc + V = {] on dVl. 

Assuming a solution exists, let us show that -Og + u is positive in 51 \ E. On the one hand, for x 
near E^, there exists some i? > such that if p[x) < Rei then 

_ 2 _ 2 

Cip{x) p-1 < Ue(x) < C2p{x) P-1 (30) 

On the other hand, since v e C^'"(il \ E), we have v{x) < cp{xY . But since v > —2/{p — 1), it 
follows that ui + V > near each E^ ; by the maximum principle we see that ui + v > in all of 
f2 \ E, and hence Ui + w is a positive solution of ([^) which is singular at all points of E. 

It remains to prove the existence of a solution to (Eq). We shall first treat the case where E 



2p 



is finite. First observe that ||/e||o,a,i^-2 < Cep ^ by Lemma Let A denote a common upper 
bound for the norm of Gg for eo sufficiently small. Then 

||G,Ml2,a,. < ACeo (31) 
In view of this estimate, we shall work in the ball 

Bieo,P) = {veCl-"{n\^):\\v\y^o.,.<Pe, -^}. 

We have already shown that Ge(/c) G B{£o,(3), for /3 large enough. At this point we shall fix 
(3 large enough so that Ggf^ G B{eQ, (3/2). We shall now make a further restriction on the e^, 
namely that, for some fixed constant a G (0, 1), we have 

aeo < ei < eo (32) 
for i = 1, . . . , K . This new restriction appears to be needed in the proof of the following Lemma. 



Lemma 9 There exists some constant c > independent o/ eo ^ 1 such that 

WQM ~ Q{vi)\\o,a,i^-2 < ^ 11^2 - Vl\\2,a,^, 

for all vi,V2 G B{eo,f3). In particular, taking V2 — we see that GgQ{v) G B{eo, (3 /2) and hence 
that the operator K. defined by IC{v) — —Gi{Q{v) + fg) maps B{eo, (3) to itself, and is a contraction 
on this hall for eo sufficiently small. 

Proof : We first establish that there exists some r > 0, independent of eo ^ 1, such that for any 
V G S(eo,/3) we have 

X G U^iB(.T,,r) =^ \v{x)\ < -ug{x). 
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Indeed, by Lemma |l|, we know that there exist constants Ci, C2 > and a radius R > such that 
cijxpp^ < Ug.{x) < C2\x\^'p^ if \x\ < Rei, 

The claim follows at once from these estimates and the fact that v G B(eo,/3), so that 

\vix)\ < pe^-^p'^ix) < cPe^-^p'^ix) 

by (H). We also note that 
for all X € ft. 

Since I^/ujI < 1/4 in each B(xi,T), we may use a Taylor expansion to obtain 

\Q{V2) - Q{V1)\{X) < C\U,\P-^{X)\V2 - 

So for X G B{xi, t) we have 

Pixr-IQM - QM\{X) < Cp(x)™f-^+^'^||«2 - ^^lIlL,. 

2 2 A''— ^ 

= Cp{x)''+—\\V2 - Willie,, < Cr'^+ — /3eo \\V2 - Wl||2,a,i.. 

The coefficient here may be taken as small as desired by choosing eo sufficiently small. Outside 
the union of these balls we use the estimates 

< CEg and also \v{x)\ < ce^ "'^ , 

where the constant c > depends on t but not on v. For p{x) > t we can neglect all factors 
involving p{x), hence for all x G J7r we have 

p{xf-''\Q{v2) - QMlix) < C{ur' + \vr')\v2 - Vi\ 

, ^ (p-mN-^) P {N-2)p-N.. II 

< Ceo h2-wi|<— e5 11^2 - wi||2,Q,i- 

for any constant B > 0, provided eg is chosen small enough. 

We also need to estimate the Holder norm of Q{vi) — Q{v2). By the fourth part of Lemma ||, 
it suffices to estimate p^^'^\V{Q{vi) — (5(f2))|. We shall only sketch this briefly. First, in each 
B{xi,T) we compute that 

VQ{v) - p + v)P-^ - wf-i -{p- l)uf-^) (VS.) + p ((u,- + v)P-^ - ur')(Vi')) . 

For X E fir there is a similar expression, except that we must be more careful about the absolute 
value signs. However, in we estimate each term individually and conclude that 



p3--|VQ(«i) - VQ(i;2)| < Cei''-^^''-''-'\\v2 - «i||2 



a, I' 



there. As before, this provides the proper estimate in this set. In B{xi,T) we use a Taylor 
expansion in the expression above to conclude the proper estimate. This completes the proof of 
the Lemma when E is discrete. 

To carry out the proof in the more general case, only fairly minor changes need to be made 
in the argument above. The most important one, already observed in §4, is that the weight 
parameter v must now be constrained to lie in the smaller interval 

_L < ^ < nrin{— ^ + 1 ^ (33) 
p— I 1 P— 1 
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Furthermore, the estimate in Lemma ^ is shghtly weaker than that in Lemma hence we must 
replace the exponent N — in the above argument by 



p-3 

q = - j^. 

p~l 

Thus we show that /C is a contraction of the balls 

The details of the argument need very little change since the nonlinear term Q{v) is the same as 
before. This completes the proof in all cases. □ 

Using this Lemma and the remarks preceding it, we now see that there exists a unique solution 
V of the equation 

V = -G, (^{\v + u,\P - u?-pul-^w) + /,) 

in B{eo, /?) for every cq small enough, and for every if-tuple e satisfying (|3^). Fixing some eo which 
is small enough for this argument to work, and replacing eo by e < eo in this whole argument 
shows that we can find a solution which is a small perturbation of Ue for every e and Ug for which 
ae < ei < e < eo. 



10 The proofs completed 

We have now completed the proof of the existence of solutions for Theorems ^ and ^. In this 
section we indicate the arguments leading to the deformation space of solutions for either problem. 
We also indicate the small changes needed to complete the proof of Theorem 3. 

10.1 The deformation spaces 

Included in the statements of Theorems |5 and ^ are the assertions that the space of all solutions 
to the equation on a domain C M with isolated singularities at {xi, . . . ,xk} C is if 
dimensional, and in general that if any component Si of the singular set is positive dimensional, 
then there is an infinite dimensional family of solutions with the same singular set. We recover 
both of these statements using the implicit function theorem. The latter case is analogous to the 
one studied in [|| , and the former mirrors the situation in . 

In both cases, the main point is that if u is the positive solution to (|l|) constructed by the 
procedure of this paper, then the linearization L — A + pu^^^ to the equation at u is nondegen- 
erate in the sense that it is surjective as a map from C^'!^{Q, \ S) to C'^'^2(^ \ This means, 
using the implicit function theorem, that all solutions to the equation which can be obtained as 
perturbations off of u by terms growing no faster than p'^ are parametrized locally by elements of 
the nuUspace of this linearization in p. We shall prove this surjectivity in the next proposition, 
or rather, we will prove the dual statement as in §6 that the linearization is injective on C^^/p. 

Proposition 7 If u ~ + v is a solution of as constructed in the last section, with all 
ti £ (ae, e), and v G C^p(r2 \ S) fl B{e,l3), then if e is sufficiently small, the linearization — 

A + puP~^ is injective as a map on C^^''^{Q, \ E). 

Proof: Suppose not, i.e. suppose there exists some sequence e^^^ with e'f'' G (ae*-^-*, e*-^') for some 
sequence e^^' 0, and a solution = u^(i) + v such that the corresponding linearization Lu^^^ 

is not injective. Thus there exists some 0^^^ e ^'f.v such that Lu^^-^cj)'^^'' = 0. We make exactly 
the same rescalings as in the proof of Theorem pi to arrive at the same contradiction. (Note, 
however, that since the el^"* must all tend to zero simultaneously, the proof simplifies somewhat.) 
The only point that needs checking is that the operator obtained in the limit of these rescalings 
is just Li or Li, according as whether S is discrete or not. The rescaling sending Ue{r) to ui(r) 
is ui{r) = e^/'^P~^^Ue{er). Applying this same rescaling to v{r) (suitably translated so that the 
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origin corresponds to the appropriate point of S) yields the function i'c(r) — e^/^P ^^v[er). Since 
\v(r)\ < Pe'^p'^ for some positive exponent q, we see that |f£(f)| < f3e'''^'^'^'p^ p'^ in some ball -Bo-/€- 
Since i/ + > and q > 0, this term tends to zero with e. Hence the limiting operator is just 

A +pu^^^, as desired. As noted above, this is a contradiction, so the proposition is proved. □ 

The next step is to compute the size of the nuUspace of Lu on Cl '^iVL \ S). We have already 
indicated that when some is positive dimensional this nuUspace is infinite dimensional, and 
hence there is an infinite dimensional space of solutions to the nonlinear problem. But when E 
is a finite set, then L is Fredholm and the nuUspace is finite dimensional. We may calculate its 
dimension using a relative index theorem. 

Proposition 8 If is discrete, this nuUspace of L„ is K dimensional. If any has positive 
dimension, then the nuUspace is infinite dimensional. 

Proof: The latter statement, that the nuUspace is infinite dimensional when some dimE; > 
follows from the theory of In fact, this nuUspace is parametrized by functions (of a certain 
negative distributional order!) along each of positive dimensional components of S. This infinite 
dimensionality can also be seen as a consequence of the fact that the globalized linearization Li 
from §6.4 has infinite dimensional nuUspace, as was pointed out in Proposition |^. 

That the nuUspace of L„ on C^'^ is K dimensional when S is discrete is proved in exactly 
the same manner as the analogous fact in 0|. We sketch this argument briefly here. Once again, 

techniques are used. Let ind (J) denote the index of L„ on p^L^. From the self-adjointness 
of La on L^(ri \ E) (with respect to Euclidean measure), we find that ind(— (5) = — ind((5), for 
any S ^ {^f}- Although the index itself requires global information to calculate, it is well known 
that the 'relative index,' i.e. the difference ind ((5') — md{d") depends only on asymptotic data 
at the points of S. Using the relative index theorem proved by Melrose, we see that the relative 
index ind {—S) — ind (S) equals 2K (the details are written out carefully in |^]. Combining these 
two facts, we find that ind (— <5) = K . However, since we have proved that L„ has no cokernel as 
a map from C^^, we conclude that the nuUspace of L„ is K dimensional. □ 

Using these facts, and a standard implicit function theorem argument as in [Q, we conclude 
the following result: 

Theorem 7 When T, is a discrete set, then the solution space to the equation ^ is locally a K 
dimensional real analytic variety. The solutions constructed in this paper lie in the smooth set of 
this variety. 

The fact that this solution space is a real analytic variety (possibly with singularities) may be 
deduced using the Ljapunov- Schmidt reduction argument, as in Q. But since we have produced 
smooth points in this variety, having found solutions u for which the corresponding linearization 
Lu is surjective, we conclude that the top stratum of this variety is K dimensional, hence that 
almost every solution is nondegenerate in this sense. 

Finally, we may sharpen the deformation result when E is not necessarily discrete, applying 
the implicit function argument from to conclude 

Theorem 8 Let <Z be any union of C'^'" submanifolds E^ of dimensions ki satisfying the 
restrictions of Theorem ^. Then the equation ^ has an infinite dimensional family of solutions. 



10.2 The singular Yamabe problem on manifolds of positive scalar cur- 
vature 

The modifications in the arguments of this paper required to solve the singular Yamabe problem 
on an arbitary compact manifold (M,go), where R{go) > 0, with singular set prescribed on an 
arbitrary finite disjoint union of smooth submanifolds of dimensions greater than zero and less 
than (or equal to) {n — 2)/2 are very minor. The equation that now must be solved is 

n — 2 7x — 2 Ti+2 
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z; > on M\T., sing(w) = I]. 

The operator in the first two terms here, i.e. the hnear part, is cahed the conformal Laplacian Lg^ 
associated to the metric qq. The Unearization of this operator around the approximate solution 
Ug is 

4(n — 1) 

Now, the rescalings of this hnearization may be effected in exactly the same ways in local 
Fermi coordinate systems around the submanifolds S^. The extra term, of order zero, in this 
operator disappears in these rescalings; the resulting model operator we need to study is exactly 
the operator Li from §6. Thus we prove, as before, that this linearization is uniformly surjective, 
provided all components of e are small enough. The fixed point argument shows that (|2^) has 
a solution v. The only point of the whole argument that needs special comment is to indicate 
where we use the assumption that the conformal class of the metric go is positive. This is when 
we show that + v remains positive on all of M \ E; at this step we use that A — R{9o) 
satisfies the maximum principle - and this is only true if R{go) is nonnegative. This completes 
the proof of Theorem |[ 
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